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Abstract 

The six-vertex model with domain wall boundary conditions, or square ice, is considered for 
particular values of its parameters, corresponding to 1-, 2-, and 3-enumerations of alternating 
sign matrices (ASMs). Using Hankel determinant representations for the partition function and 
the boundary correlator of homogeneous square ice, it is shown how the ordinary and refined 
enumerations can be derived in a very simple and straightforward way. The derivation is based 
on the standard relationship between Hankel determinants and orthogonal polynomials. For the 
particular sets of parameters corresponding to 1-, 2-, and 3-enumerations of ASMs, the Hankel 
determinant can be naturally related to Continuous Hahn, Meixner-PoUaczek, and Continuous 
Dual Hahn polynomials, respectively. This observation allows for a unified and simplified 
treatment of ASMs enumerations. In particular, along the lines of the proposed approach, 
we provide a complete solution to the long standing problem of the refined 3-enumeration of 
AMSs. 
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1. Introduction 



The six-vertex model on a square lattice with domain wall boundary conditions (DWBC) was 
introduced in [1] and subsequently solved in [2], where a determinant formula for the partition 
function was obtained and proven (see also [3]). Analogous determinant formulae has been 
given also for the boundary one point correlation functions [4]. The model, in its inhomoge- 
neous formulation, i.e., with position-dependent Boltzmann weights, was originally proposed 
within the theory of correlation functions of quantum integrable models, in the framework of 
the quantum inverse scattering method [5]. The model was later found to be deeply related 
with the problems of enumeration of alternating sign matrices (ASMs) [6-11] and domino 
tilings [12-15]. It should be mentioned that ASM enumerations appear to be in turn deeply 
related with quantum spin chains and some loop models, via Razumov-Stroganov conjec- 
ture [16]; for recent works, see for instance Ref. [17] and references therein. 

In its homogeneous version, the six-vertex model with DWBC admits usual interpretation 
as a model of statistical mechanics with fixed boundary conditions, and it may be seen as a 
variation of the original model with periodic boundary conditions [18-21]. The latter was 
originally proposed as a model for two-dimensional ice (hence the alternative denomination: 
'square ice'), and has been for decades a paradigmatic one in statistical mechanics [22]. Till 
now, specific results for the six-vertex model with DWBC at particular values of its parameters 
(in fact, mainly derived within the context of ASMs) were obtained from general results for 
the inhomogeneous version, first specializing the parameters to the considered case, and then 
performing the homogeneous limit at the very end, and hence once for each particular case. 
Each time, the homogeneous limit was an hard task on its own right, and a specific approach 
was devised to work it out in each single case [10, 23, 24]. 

The purpose of the present paper is to explain how one can proceed the other way around, 
first performing the homogeneous limit once for all for the model with generic vertex weights, 
and then specializing the result to the case of interest. The homogeneous limit has already been 
done in [3] for the partition function, and in [4] for the boundary one point correlation func- 
tions, resulting in Hankel determinant representations for these quantities. Here we show how, 
by specializing the parameters to some particular values, and exploiting the standard relation- 
ship between Hankel determinants and orthogonal polynomials, all previously known results 
concerning ordinary and refined enumerations can be derived in a very simple and straight- 
forward way. It appears that for the particular sets of parameters corresponding to 1-, 2-, and 
3-enumerations of ASMs, the Hankel determinant is naturally related to Continuous Hahn, 
Meixner-PoUaczek, and Continuous Dual Hahn polynomials, respectively. The approach which 
we propose here allows for a unified and simplified treatment of ASMs enumerations, including 
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the problem of refined 3 -enumeration. 

The paper is organized as follows. In the next Section we recall the definition of the model 
and some related results, with particular attention to the Hankel determinant representations 
for the partition function and the one-point boundary correlator. In Section 3, known results on 
ASMs are reviewed, together with their connection with the square ice. In Section 4, we show 
how the Hankel determinant entering the partition function of square ice can be reinterpreted, 
in the three cases corresponding to 1-, 2- and 3-enumeration of ASMs, as the Gram determinant 
for an appropriate choice of Continuous Hahn, Meixner-PoUaczek, and Continuous Dual Hahn 
polynomials, respectively, thus providing a very simple and straightforward derivation of the 
known results for ASM enumerations. In Section 5 we exploit the fact that these polynomials, 
being of hypergeometric type, satisfy some finite difference equations, which can be translated 
into recurrence relations for the one-point boundary correlator. These recurrences in turn can be 
solved, thus giving the results for the refined 1-, 2- and, especially interesting, 3-enumeration 
of ASMs. We conclude in Section 6 with a discussion of the proposed approach. 

2. Square ice with DWBC 

Let us start with recalling the formulation of the model. The six-vertex model, which was 
originally proposed as a model of two-dimensional ice, is formulated on a square lattice with 
arrows lying on edges, and obeying the so called 'ice-rule', namely, the only admitted configu- 
rations are such that there are always two arrows pointing away from, and two arrows pointing 
into, each lattice vertex. An equivalent and graphically simpler description of the configura- 
tions of the model can be given in terms of lines flowing through the vertices: for each arrow 
pointing downward or to the left, draw a thick line on the corresponding link. The six possi- 
ble vertex states and the Boltzmann weights wi assigned to each vertex according to its state 
i {i = 1, . . . , 6) are shown in Fig. [T] As it has been already stressed in the Introduction only 
the homogeneous version of the model, where the Boltzmann weights are site independent, is 

+ + + + + + 
+ + + + + + 

Figure 1: The six allowed types of vertices in terms of arrows (first row), in terms of lines 
(second row), and their Boltzmann weights (third row). 



3 



(a) 



(b) 



Figure 2: One of the possible configurations of the model with DWBC, in the caseA^ 
in terms of arrows; (b) in terms of lines. 



5: (a) 



considered here. 

The DWBC are imposed on the N x N square lattice by fixing the direction of all arrows 
on the boundaries in a specific way. Namely, the vertical arrows on the top and bottom of 
the lattice point inward, while the horizontal arrows on the left and right sides point outward. 
Equivalently, a generic configuration of the model with DWBC can be depicted by N lines 
flowing from the upper boundary to the left one. This line picture (besides taking into account 
the 'ice rule' in an automated way) is intuitively closer to ASMs recalled in the next Section. 
A possible state of the model both in terms of arrows and of lines is shown in Fig. 2. 

The partition function is defined, as usual, as a sum over all possible arrow configurations, 
compatible with the imposed DWBC, each configuration being assigned its Boltzmann weight, 
given as the product of all the corresponding vertex weights. 



Zn — 



(2.1) 



arrow configurations i=l 
with DWBC 



Here denotes the number of vertices in the state i in each arrow configuration (ni H hng = 

The six-vertex model with DWBC can be considered, with no loss of generality, with its 
weights invariant under the simultaneous reversal of all arrows. 



wi — W2 —'■ a, ws — W4 —: b, — wq —: c. 



(2.2) 



Under different choices of Boltzmann weights the six- vertex model exhibits different behaviors, 
according to the value of the parameter A defined as 

+ 6^ - 



2ab 



(2.3) 



It is well known that there are three physical regions or phases for the six-vertex model: the 
ferroelectric phase, A > 1; the anti-ferroelectric phase, A < — 1; and, the disordered phase. 
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— 1 < A < 1. In the present paper we shall discuss only some particular cases, with values of 
the Boltzmann weights that correspond to the disordered phase. A convenient parametrization 
of the Boltzmann weights in this phase is 

a = sm{\ + ri), 6 = sin(A — 77), c = sin2?7. (2.4) 

With this choice one has A = cos 2r]. The parameter A is the so-called spectral parameter and 
1] is the crossing parameter. The physical requirement of positive Boltzmann weights, in the 
disordered regime, restricts the values of the crossing and spectral parameters toO<r7<7r/2 
and T] < \ < 7r/2 — rj. 

An exact representation for the partition function was obtained in Ref. [3]. When the 
weights are parameterized according to (12.41) such representation reads 

Z, = MA-,)sin(A + ,r- ^ ^2.5) 

n ik\y 

k=l 

where Z is an x matrix with entries 

2;.^ = .^ si^^^ (2 6) 

^ 9A-?'+^ sin(A — 77) sin(A + rj) ' 

Here and in the following we use the convention that indices of N x N matrices run over the 
values j, A; = 0, 1, . . . , A^ — 1. 

This formula for the partition function has been obtained as the homogeneous limit of a 
more general formula for a partially inhomogeneous six-vertex model with DWBC. The inho- 
mogeneous model, with site-dependent weights, is defined by introducing two sets of spectral 
parameters {\j}jLi and {z/^}^^, such that the weights of the vertex lying at the intersection of 
the j-th column with the fc-th row depend on \j — rather than simply on A, still through for- 
mulae (12.41) . The inhomogeneous model, though apparently more complicate, can be fruitfully 
investigated through the Quantum Inverse Scattering Method, see papers [1-3] and book [5] 
for details. As a result, the partition function of the inhomogeneous model is represented in 
terms of certain determinant formula which, however, requires great caution in the study of its 
homogeneous limit, Uk ^ and Xj A, since in this limit the determinant possess A^^ — A^ 
zeros that are cancelled by the same number of singularities coming from the pre-factor. A 
recipe for taking such a limit was explained in detail in Ref. [3] where formula (12.51) was orig- 
inally obtained. Subsequently, formula (12.51) was used in papers [25, 26] to investigate the 
thermodynamic limit, A^ ^ 00, of the partition function. In these studies the Hankel nature 
of the determinant appearing in (12.51) . a natural outcome of the homogeneous limit procedure, 
was exploited through its relation with the Toda chain differential equation or with the random 
matrix partition function. 
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The aim of this paper is to explain how this Hankel determinant formula (and its analogue 
for the boundary correlator, given below) can also be used in application to some well-known 
problems of combinatorics, such as enumerations (and the so-called 'refined' enumerations) 
of ASMs. These problems and the known results are reviewed in the next Section. It is to be 
mentioned that though in these combinatorial problems one deals in fact with the homogeneous 
six-vertex model with DWBC, the Hankel determinant formula for the partition function was 
never discussed previously in this context. Instead, the more complicated determinant formula 
for the inhomogeneous model partition function was used in the combinatorial proofs. In these 
proofs, the homogeneous limit can be regarded, in fact, as the most complicated part on the 
way to the result. Consequently, one can expect that extracting relevant information directly 
from the homogeneous model should be technically much simpler. We will show that this is 
indeed the case since some standard classical orthogonal polynomials can be naturally related 
(for some particular choices of the parameters) to the Hankel determinant in (12.51) . 

In addition to the partition function, we shall discuss here also boundary one point corre- 
lation functions. In general, two kinds of one point correlation functions can be considered 
in the six- vertex model: the first one ('polarization') is the probability to find an arrow on a 
given edge in a particular state, while the second one is the probability to find a given vertex in 
some state i. If one restricts to edges or vertices adjacent to the boundary, then such correlators 
are called boundary correlators. Following the notations of paper [4], where these boundary 

(r) 

correlators were studied, let G\f denote the probability that an arrow on the last column and 
between the r-th and the r + 1-th rows (enumerated from the bottom) points upward (or, in the 

(r) 

line language, that there is no thick line on this edge), and let Hj^' denote the probability that 
the r-th vertex of the last column is in the state z = 5 (or that the thick line flows from the top 
to the left), see Figs. [T] and |2l The first correlator, G^^\ is, in fact, the boundary polarization, 

(r) 

whose interpretation is more direct from a physical point of view, while the second one, Hj^ , 
is closely related to the refined enumerations of ASMs. It is easy to see that, due to DWBC, the 
two correlators are related to each other as follows 



In Ref. [4] both correlators were computed using Quantum Inverse Scattering Method for 
the inhomogeneous six-vertex model. In the homogeneous limit, which is the situation we are 
interested in here, determinant formulae generalizing (12.51) were found for these correlators. 

(r) 

For instance, for H]^ , the following expression was derived 




(2.7) 




(iV- 1)! sin277 



det n 



(2.8) 



[sin(A — 77)] ^ [sin(A + ■r])~\ 



^-^+1 det Z 
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where the N x N matrix H differs from Z only in the last column, 



n 



jk 



[sine]'' \mi{e — 2r] 



N-r 



[sm{e + X — T] 



for k = 0,...,N -2 
for A; = - 1 



(2.9) 



e=0 



A similar expression is valid for G^^^ as well. In what follows we shall focus on H^i!f'; the 
results for Gj^-* will follow immediately from relation (I2.7I) . From the DWBC it immediately 



(r). 



follows that G^^^ = 1, and therefore, from (I2.7I) . correlator H'^^ has to satisfy 



(r) 

N 



N 



r=l 



(r) 
N 



1. 



(2.10) 



In what follows this normalization condition will be used in application to the generating func- 
tion ofH^. 



3. Alternating Sign Matrices 



An alternating sign matrix (ASM) is a matrix of I's, O's and — I's such that in each row and 
in each column ( i) all nonzero entries alternate in sign, and ( ii) the first and the last nonzero 
entries are 1. An example of such matrix is 

/O 1 o\ 

1-110 

1 0. (3.1) 

10 
\0 ly 

There are many nice results concerning ASMs, for a review, see book [11]. Many of these 
results have been first formulated as conjectures which were subsequently proved by different 
methods. 

The most celebrated result concerns the total number A(N) of N x N ASMs. It was 
conjectured in papers [6,7] and proved in papers [9, 10] that 

^ (3fc-2)!(A;- 1)! ^ 



A{N) = n 



k=l 



{3k-2)\ 

{2k - ly. {2k - 2)\ {2N-k)\ 



n 



(3.2) 



Other results concern the weighted enumerations or the so-called x-enumerations of ASMs. 
In x-enumeration the matrices are counted with a weight x'' where k is the total number of 
— 1 entries in a matrix (the number x here should not be confused with the variable x widely 
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used in the following Sections). The number of x-enumerated ASMs is denoted traditionally 
as A{N; x). The extension of the x = I result above to the case of generic x is not known, but 
for a few nontrivial cases, namely x = 2 and x = 3, closed expressions for x-enumerations are 
known (the case x = is trivial since assigning a vanishing weight to each —1 entry restricts 
the enumeration to the sole permutation matrices: A(n] 0) = n\). The result in the case x = 2, 
related to the domino tilings of Aztec diamond [12, 13, 15], have been obtained in [7], and reads 

A(iV; 2) = 2^(^-1)/^ (3.3) 

The answer in the case of 3-enumeration, again conjectured in [6,7], was subsequently [10] 
proved to be 

yl(2m+2; 3) = 3^ (3^ + 2)!^! ^^^^^^ 
[(2m+l)!]' 

(3.4) 

Another class of results concerns the so-called refined enumerations of ASMs. In the refined 
enumeration one counts the number of x ASMs with their sole 1 of the last column at 
the r-th entry. The refined enumeration can be naturally extended to be also an x-enumeration. 
The standard notation for the refined x-enumeration is A(N, r; x); in the case x = 1 one writes 
simply A{N, r) just like A{N) for the total number of ASMs. It was again conjectured in [6,7], 
and proved in [23] that refined enumeration of ASMs is given by 

AiN,r)= ^^-]Jl^^~' ^ AiN). (3.5) 

In the case of the refined 2-enumeration it has been shown in [7, 12, 13, 15] that 

^(^>^;2) = ^(^~^^)A(Ar;2). (3.6) 

The case of the refined 3-enumeration appears to be much more complicate. Direct computer 
enumeration does not suggest any factorized form and no conjecture concerning A{N, r; 3) has 
ever been proposed. Recently, Stroganov obtained a certain representation for the correspond- 
ing generating function [27], allowing in principle a recursive computation of the numbers 

A{N,r-3). 

The most direct way to recover all the previous results is based on a nice bijective corre- 
spondence between ASMs and six-vertex model, which has been pointed out in [8, 12, 13], and 
applied for the first time in [10], namely, that to each single N x N ASM corresponds one and 
only one arrow configuration of the six- vertex model on the N x N square lattice with DWBC. 
The correspondence between matrix entries and vertices is depicted in Fig.|3l As an example. 



A(2m+1;3) = S^^^+i) JJ 



7 1 



(3A;- 1)! 
{m + k)\ 
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+ + + + + + 

1 -1 

Figure 3: Vertex states — ASM's entries correspondence. 



matrix (13.11) corresponds to the configuration of Fig.|2land vice versa. 

As an immediate consequence of this correspondence, ASM enumeration is exactly given 
by the partition function of square ice, when all vertex weights are set equal to unity. More 
generally, the number of — I's in a given ASM being equal to the number of vertices of type 
6 (see Fig. |3l), and the number of vertex of type 5 and 6 being constrained by the condition 

— uq = N, we readily get 



A(iV;a;) = (l-x/4)-^'/2^-^/2^ 



N 



X=n/2 

T;=arcsin(Y^/2) 



(3.7) 



Therefore, x-enumeration of ASM corresponds to the computation of the partition function 
of square ice on the subset of parameters space given by a = 6. In this correspondence, 
values of x belonging to the interval (0, 4) corresponds to the disordered regime of the model, 

-1 < A < 1. 

This nice correspondence can be further extended to the refined x-enumeration of ASMs. 
In the language of square ice, the ratio A{N, r; x) /A(N; x) can be rephrased as the probability 
of finding the unique vertex of type 5 on the right boundary at the r-th site, which is exactly the 

(r) 

definition of the boundary correlator H}^' . Explicitly, one has 



A(N,r;x) 
A{N;x) 



H 



(r) 
TV 



A=V2 • (3-8) 

»7=arcsin ( •y/2:/2) 



Thus, being able to compute the partition function and the boundary correlator for some partic- 
ular choice of parameters, one immediately obtains x-enumerations and refined x-enumerations 
of ASMs respectively, for some corresponding values of x. 



4. The partition function and enumerations of ASMs 
4.1. Preliminaries 

From representation (12. 5I) . it is evident that the evaluation of the the partition function in the 
homogeneous case essentially reduces, modulo a trivial pre-factor, to the calculation of the 
determinant of the N x N matrix Z, which is a Hankel matrix, whose entries do not depend 
on A^. There is a standard method to treat such determinants, which is based on the theory of 
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orthogonal polynomials [28, 29], and has proven to be quite powerful when some assumption 
are verified. 

Let us assume that the entries of our Hankel matrix are written in the canonical form 



oo 



Zjk= / x^^^^i{x)dx. (4.1) 



oo 



Let us moreover suppose that there exist a complete set of polynomials {Pn{x)}^^Q subject to 
the orthogonality condition 



Pj(a^) Pkix) j^{x) dx = hjdjk ■ (4.2) 
Then, denoting by k„ the leading coefficient of Pri(x), 

Pn{x) = + . . . , K„ 7^ , (4.3) 

and using standard properties of determinants and of orthogonal polynomials, one obtains for 
the determinant of the N x N matrix Z the following formula 



detZ = det 

7V-1 



Pj{x)pk{x)iJ,{x)dx 



(4.4) 

j,k=0 



n=0 " 

Obviously, this formula may turn out useful provided that the set of polynomials which are 
orthogonal with respect to the weight fi{x) can be identified. 

In the case of matrix Z with entries (12.61) one can easily obtain the weight fi{x) using the 
representation 

= r e(^-/2)- ^^^^^^ ax (4.6) 

sin(A — rj) sin(A + rj) sinh fx 



which is valid for A and r] corresponding to the disordered regime [26]. Unfortunately, appro- 
priate polynomials are not available in general and the previous scheme cannot be fulfilled for 
generic values of A and r]. However, for some very particular values of these parameters the ap- 
propriate orthogonal polynomials appear to be known, and have just to be suitably chosen in the 
framework of Askey scheme of hypergeometric orthogonal polynomials [30]. To be precise, 
there are essentially three cases, indicated in the phase diagram of the model, see Fig.|4l which 
fit into the scheme: i) the so-called 'free fermion' line 77 = 7r/4 and n/A < A < 37r/4; ii) the 
A = 1/2 symmetric point (or 'ice point') r] = n/6 and A = 7r/2; and Hi) the A = —1/2 sym- 
metric point 1] = 7r/3 and A = tt/2. With 'symmetric' here we mean that these points lie on the 
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h/c 




1 







a/c 



Figure 4: The phase diagram of the model, with ferroelectric (F), antiferroelectric (AF) and 
disordered (D) phases, separated by the solid lines. The three considered cases, all belonging 
to the disordered phase, are shown in bold: the free fermion line, and the three points corre- 
sponding to 1-, 2- and 3-enumeration of ASMs. 

line a = b. The polynomials corresponding to these three cases are Meixner-PoUaczek polyno- 
mials. Continuous Hahn polynomials and Continuous Dual Hahn polynomials, respectively. In 
the rest of this Section we give details of computation for each of these cases. 

It is worth mentioning that all previously listed choices of parameters, for which the Hankel 
determinant appearing in (12.51) turns out to be related to some set of classical orthogonal poly- 
nomials, exactly cover 1-, 2- and 3-enumerations of ASMs. On the other hand, the fact that no 
set of polynomials in the Askey scheme corresponds to some other choice of the parameter r] 
(even with A set equal to 7r/2) is likely to be deeply related to the lack of factorizable formulae 
for x-enumerations of ASMs other than for x = 1,2, 3. 

4.2. The free fermion line 

We shall start with the r] = tx /A case which is technically the simplest one, though the parameter 
A is not fixed, 7r/4 < A < 37r/4. In this case the orthogonality weight ^{x) can be written as 




[30] 




where 
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we find that the polynomials related to the Hankel determinant in this case are 

P.(x) = Pi^/^)(|;2A-|). 

Using, for convenience, the parameter = 2A — 7r/2 G (0, tt) instead of A, we have 

2 



(4.10) 



sm( 



2"n! 



(4.11) 



Inserting these expressions into (14.51) . we straightforwardly obtain the following result for the 
partition function 



smffl -n- 2 



n 



2'^n! 



sm( 



1. 



Using Eqn. (I3.7I) . we immediately recover formula (13.31) for 2-enumeration of ASMs. 



(4.12) 



4.3. The ice point and the total number of ASMs 

At the ice point, or A = 1/2 symmetric point, the computation is very similar. For r] = ir/G 
and A = 7r/2 the orthogonality weight reads 



/i(x) 



sinh fx 1 

o 

sinh Air'^ 



^, 1 x\^f2 X 
V[- + i- r - + 
'3 6 / V 3 6 



(4.13) 



where the triplication formula for the F-function 



33X-1/2 / 1\ / 2 

r(3x) = T{x)t[x + - j Ff X + - 



(4.14) 



has been used. The orthogonality condition for Continuous Hahn polynomials is 

— / Pn{x'-, a, b, c, d)pm{x; a, b, c, d) F(a + ix) T(b + ix) F(c — ix) T{d — ix) dx 

F(n + a + c) T{n + a + d) V{n + b + c) V{n + b + d) 



{2n + a + b + c + d - l)V{n + a + b + c + d - I) n\ 
and the polynomials are given by 

/ , „ .^{a + c)n{a + d)n ^ f -n, n + a + b + c + d-l, a + ix 

Pn{x;a,b,c,d) =1 ■ 3i^2 , , , 

n\ \ a + c , a + d 



Km (4.15) 



1 . (4.16) 



Orthogonality condition (14.151) is valid if the parameters a, b, c, d satisfy Re (a, b, c, d) > 0, 
a = c and b = d. Comparing (14.131) with (14.151) we are naturally led to the choice of parameters 
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a = c = 1/3 and h = d = 2/?>. Hence, the appropriate polynomials to be associated to the 
Hankel determinant in this case are 



'x 1 2 1 2 



i"(2/3)„ 3i^2 



— n, n + 1, l/3 + ia;/6 



2/3, 1 

For the normalization constant and the leading coefficient we have the expressions 

2(3n+l)! (2n)\ 



(4.17) 



(4.18) 



" (2n + 1)33"+V2n!' 6"(n!)2 
where the triplication formula, Eqn. (14.141) has been used to simplify the normalization constant 

hn. 

Substituting the obtained values of hn and /t„ in expression (14.51) for the determinant, taking 
into account the value of the prefactor in (12.51) for A = tx /2 and rj = n/6, and cancelling 
whatever possible, we arrive to the following value for the ice point partition function 

2 



X=tt/2 
r)=7r/6 



n 

n=0 



(3n + l)!n! 
(2n)! (2n + 1)!' 



(4.19) 



The product expression here gives exactly the total number of ASMs, A{N), since by formula 
(13.71) the first factor relates A(N) with the partition function. 



A{N) = (3/4)-^'/2^jv 



A=7r/2 • 
??=7r/6 



(4.20) 



Thus, we have easily recovered the celebrated result, Eqn. (13.21) . for ASMs enumeration directly 
from the Hankel determinant formula (12.51) . 



4.4. The A = —1/2 symmetric point and 3-enumeration of ASMs 

In this case, A = 7r/2 and r] = n/S, the weight yu(x) can be rewritten, using (14.141) . in the form 

1 



/i(x) 



sinh |x 
sinh |a; Svr^ 



r(if)r(| + i|)r(| + if; 



r(i|; 



(4.21) 



This expression recalls the weight for Continuous Dual Hahn polynomials S'„(a;^; a, b, c), which 
are defined by 

' —n , a + ix , a — ix 
a + b , a + c 

and, for real and nonnegative values of parameters a, b, c, satisfy the orthogonality condition 



Sn{x'^; a, b, c) = (a + 6)„(a + c)„ 3F2 



(4.22) 



27r 



S.m{x'^; a, b, c)Sn{x'^; a, b, c) 



r(a + ix)r(6 + ix)r(c + ix) 



r(2ix) 



dx 



r(n + a + b)T{n + a + c)T{n + b + c) n\ Snm ■ (4.23) 
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At a first glance, however, Continuous Dual Hahn polynomials do not seem appropriate for the 
evaluation of the Hankel determinant: on one hand, these are polynomials in rather than in 
X, and on the other hand the integration domain in (14.231) is restricted to the positive half-axis. 
These inconvenience may nevertheless easily be circumvented, thanks to the following nice 
feature of matrix Z, holding for A = 7r/2 and generic rf. whenever the sum of indices j ' + A; is 
odd, the corresponding entry Zjk vanishes. In other words, the matrix Z at A = 7r/2 exhibit 
a chessboard pattern of vanishing and non vanishing entries. As can be seen from the Laplace 
formula for the determinant of the sum of two matrices, this chessboard structure immediately 
implies that the determinant of such matrix always factorizes onto two determinants of smaller 
matrices with no vanishing entries. 
To be precise, let us denote 



D 



N 



detZ 



A=7r/2 
r;=7r/3 



det 



AT-l 



j,k=0 



(4.24) 



where |u(x) is given by (14.211) . Since the orthogonality weight is even, ii{x) = x), one has 



x^'^'' /x(x) dx 



and consequently the following factorization arises 




fi{x) dx ifj + k is even 
ifj + k is odd 



D. 



2m 



D 



2m+l 



(4.25) 



(4.26) 



Here Dm and Dm are determinants of m x m Hankel matrices, built only from the even 
moments of the weight /i(x), 

Z^L^) = det 

where 



m—l 



a 



0,1 



(4.27) 



J i,fc=0 



^i^^\x) = 2n{x), /i(^)(x) = 2x^fi{x). 



Let pi^'' (x) be the polynomials subject to the orthogonality condition 

) 

pf\x^)pi^\x^)fi^''\x) dx = hf6,k, 



(4.28) 



(4.29) 



then, in analogy with Eqn. (14.51) . we have the following formula 



D^^^ = det 



^pf\x')pi^\x')f^^^\x)dx 



-1 m—l 



'kI 



J j,k=0 



m-1 



n 



(4.30) 



n=0 L"-" 
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where k^T^ is the leading coefficient of the polynomial pi'^''(a;^). Thus, Continuous Dual Hahn 
polynomials can be related to the Hankel determinants D^m by properly specializing the pa- 
rameters (a, 6, c) in each of the two cases, cr = and cr = 1. 

Comparing Eqns. (I4.28I) and (I4.21I) with the orthogonality weight (I4.23I) . we are led to 
specialize the parameters (a, 6, c) to the values (0, 1/3, 2/3) for a = and (1, 1/3, 2/3) for 
(T = 1. Thus, we have 



Fn \ 1 " > 36' ' 3 3 



3-^2 



-n , ia;/6 , —ix/Q 
1/3, 2/3 



and 



3-^2 



— n , l + ix/6, 1— ia;/6 
4/3 , 5/3 

For the normalization constant and the leading coefficient we have 



n\{3n)\ 

'33n-l/2 



and 



' n " 33n-l/2 • 



36 



36 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



respectively. The F-function triplication formula, Eqn. (14.141) has again been used to simplify 
the expressions for the normalization constants hn"* and hn"^ . Using (14.301) we obtain 

m— 1 



fc=0 
m— 1 



(4.35) 



fc=0 



Substituting these values into (I4.26I) . taking into account the pre-factor of Eqn. (12. 5I) . and 
cancelling whatever possible, we obtain the following expression for the partition function 

at A = 7r/2 and 77 = 7r/6, 



■'2m 



\=-k/2 



\=tt/2 
»7=7r/3 



4m^ 



yin +m 



(3m)! 



(2m+l)2 



m—1 

n 

fc=0 

m—1 



■^2m+l 

Recalling (13.71) which now reads 

A{N; 3) = 2^'- 



r^m^+2m+l/2 JJ' 



k=0 



{3k + 2)1 ' 
{m + k)\ 

{3k + 2)1 
{m + k + 1] 



(4.36) 
(4.37) 



1 

7! 



N 



A=7r/2 
7?=7r/3 



(4.38) 
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formulae (13.41) for 3-enumeration of ASMs are readily recovered. 

In conclusion, in this Section we have given explicit formulae for the partition function of 
the six-vertex model with DWBC for some particular values of parameters A and i]. Analogous, 
essentially equivalent, expressions were already known from the investigations of ASMs x- 
enumeration. It is worth emphasizing that the approach presented here allows to recover these 
results in a very simple and straightforward way. The keystone of the whole approach is the 
nice connection with known classical orthogonal polynomials that will now be further exploited 
to explore the boundary correlation function in application to the refined x-enumerations of 
ASMs. 



5. The boundary correlator and refined enumerations of ASMs 
5.1. Preliminaries 

In this Section we shall show how the knowledge of the suitable set of classical orthogonal 
polynomials associated to each of the three considered cases, can be further exploited to de- 
rive explicit answers for the one point boundary correlation function. As a consequence the 
results for the refined a;-enumerations of ASMs can be obtained; the result for the refined 3- 
enumeration of ASMs is of primary interest since it was not conjectured previously. 

In what follows, for the sake of simplicity and clarity, we shall often ignore the overall 
normalization of the boundary correlator. As already discussed at the end of Section |2l the 
proper normalization can always be restored trough the use of Eqn. (12.101) . 

The basic idea is very simple, stemming from the fact, see Refs. [28,29], that the polynomi- 
als associated to a given orthogonality weight can be in turn represented as determinants. 



Pn-i{x) = const X det>V, (5.1) 

where N x N matrix W differs from Z, defined by (14.11) . just only in the elements of the last 
column. 



Z. 



■jk 



jk 



for A; = 0, . . . , - 2, 
X' for A; = — 1 
The boundary one point correlator, see Eqn. (12.81) . therefore reads 



(5.2) 



H 



(r) 
N 



const X 



sin(A + T]) 
sin(A — T]) 



PN-l{de 



(sine) 



r-l 



sm e 



N-r 



(5.3) 



e=0 



[sin(e + A-r])]^-! 

This representation is completely general, however it can be fruitfully exploited only when the 
explicit form of the polynomial entering it is known, which happens precisely in each of the 
three cases under consideration, as explained in the previous Section. 
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Indeed, in these three cases, the polynomials of interest, being of hypergeometric type, are 
known to satisfy the following finite difference equations with respect to their variable, see, 
e.g., Ref. [30]. Denoting y{x) = P^\x; 0) for Meixner-PoUaczek polynomials one has 

Q^'^{a — ix)y{x + i) + 2i[x cos (j) — [n + a) sin (l)\y{x) — e~"^(a; + \x)y{x — i) = 0. (5.4) 

The Continuous Hahn polynomials satisfies 

B{x)y{x + i) - [B{x) + D{x) + n{n + a + b + c + d - l)]y{x) + D{x)y{x - ia;) = (5.5) 

where y{x) = Pn{x; a, 6, c, d) and 

B{x) = (c - ix){d - ix), D{x) = (a + ix){b + ix). (5.6) 

For the Continuous Dual Hahn polynomials one has 

B{x)y{x + i) - \B{x) + D{x) + n\y{x) + D{x)y{x - ix) = (5.7) 

where y{x) = Sn{x^] a, 6, c) and 

_ {a - ix){h - ix){c - ix) _ (a + ix)(& + ix)(c + ix) 

2ix(2ix-l) ' 2ix(2ix + l) • 

The approach we shall apply here to compute the boundary correlator is based on the fact that 
each of these finite-difference equations for the polynomials can be directly translated into a 
recurrence relation for the boundary correlator which in turn can be solved explicitly. 

The derivation of the recurrence relations for the boundary correlator in the three cases is 
quite similar and will be explained below in detail in each case. The general idea underlying 
the procedure to obtain the recurrence relation is based on the simple relation y{de ± i) = 
e^'^l/(9e)e^'^ which allows us to derive from each finite difference equation a relation of the 
form 

y{d,)KJ{e) = 0. (5.9) 

e=0 

Here is some linear differential operator whose form is determined by the finite difference 
equation, and f{e) is a trial function. Under special choice of this function (15.91) becomes a 
linear recurrence relation for the boundary correlator. 

To find this function we shall use also the fact that the formula (15.31) can be further rewritten 
as follows 

= const X [pN^,{d,)[g{e)f~\uj{e)Y-^] 
where uj{e) and g{e) are given by 



(5.10) 

e=0 



^^^^ sin(A + r/) sine ^^^^ sin(A - r/) sin(£: - 2//) (5 11) 

sin(A — rf) sin{e — 27]) ' s\.n{2ri) s\.n{e + A — r/) ' 
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and related to each other as 

9{e) = (5.12) 

It is clear that in order to rewrite Eqn. (15.91) as a recurrence relation for if^'' the function f{x) 
should be chosen of the form 

f{e) = [9{e)f-\uj{e)Y~\{e) (5.13) 

where T{e) is still arbitrary. Then, rewriting the operator in terms of the differential operator 

^ „ sin e sin(e — 2?]) ^ , ,s 

D, := ud^ = r\ d,. (5.14) 

sm irj 

and reexpressing all quantities in terms of the variable uj, the expression for T{e) can be easily 
chosen in such a way that Eqn. (15.91) turns into a recurrence relation for the boundary correlation 
function. The actual procedure becomes apparent after turning to the examples. 

As a last comment here it should be mentioned that the recurrence relation for H]^' can 



equivalently be viewed as an ordinary differential equation for the generating function 

N 

Hn{z) ■.= Y,h'P z'-\ Hn{1) = 1. (5.15) 

Thus the problem of solving the obtained recurrence relation can also be regarded as that of 
finding a polynomial solution to the corresponding differential equation. In all cases considered 
below, such differential equations appear to be at most of the second order. 

5.2. The free fermion line 

We start with reminding that this is the case when t] = n/A, with the spectral parameter free to 
vary within the interval 7r/4 < A < 37r/4, or, using = 2A — 7r/2 we have < < vr. The 
value A = 7r/2 (or = 7r/2) corresponds to 2-enumerated ASMs. 

The boundary correlator in this case is given by the formula (15.101) where 

/ N / , /^N / \ sin((/)/2) cose 

^ e = - cot 0/2 tans, g{e) = . , / , (5.16) 

sm(e + 0/2) 

and PN~iix) is given by the formula (14.101) . i.e., being the particular case of the Meixner- 
PoUaczek polynomial. The finite-difference equation (15.41) in this case reads 



— cos — ( — - 1 sin ( 



-e-'M^ + f )P^-i(^-4i) = 0. (5.17) 
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Using PN-i{de ± 4i) = e^'''^pAr_i((9e)e^'''^ we derive the condition 



where K;, is the first order linear differential operator 

sin e cos e 



£ = 



sm I 



de sin(2e + (/)) + AT - 1 



(5.18) 



(5.19) 



and r(e) is an arbitrary function which will be suitably chosen below to turn equation (15.181) 
into a recurrence relation for the boundary correlator. Our aim now is to explain how this 
function can be found and the recurrence relation can be obtained. 
First, using the operator 

-Df = sin e cos ed^ = ud,,. 



and changing to the variable lu using 

sin(2e + 0) {uj - l){auj + 1) 



sin( 



a := tan^(0/2), 



we find 



au'^ + 1 

= \d,{u - l){auj + 1) + (iV - l){auj^ + 1)| ,\ . 
I J aw^ + 1 

Next, let us define the operator by 
where g is given by (15.161) . Using 

g-\D,g) : 

we find 



' UJ - 1 



D,{auo^l) -{N -I) 



uo-l 



(5.20) 

(5.21) 
(5.22) 

(5.23) 

(5.24) 
(5.25) 



auj"^ + 1 

Now the choice of the function T{e) is evident, since in terms of the operator the relation 
(15.181) reads 

PN~i{de)g''-'K,uj'-'T =0 (5.26) 

e=0 

and therefore if we choose T{e) to cancel the factor standing outside the braces in (15.251) then 
we immediately obtain 



PN^im (7^-1 [D,{auj + 1) - (iV - l)]uj''^ 



0. 



e=0 



(5.27) 



Finally, reminding that = Lod^ this last equation directly leads to the recurrence relation 



ar - (iV - r)//);^ = 



(r) 



(5.28) 
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where we have used (I5.10I) . Note, that actual choice of the function r(e) have been governed to 
have all coefficients in the braces in (15.251) to be polynomials in uj. Note also the 'anti-normal' 
ordering of the differential operator there. Below we shall proceed in other cases exactly in the 
same way. 

Solving the recurrence (15.281) we find 



" > >■ - 1 y |i + tari2(^,/2)|"-' 



Here the proper normalization is achieved by directly satisfying condition (12.101) . 

Thus, the result of Ref. [4] is readily recovered. Moreover, recalling (13.81) . and specializing 
A = 7r/2, that is = 7r/2 or tan(</)/2) = 1, the expression (13.61) for the refined 2-enumeration, 
obtained in [7, 12, 13, 15], is immediately reproduced. 

5.3. The ice point and the refined enumeration of ASMs 

Now we turn to the interesting case of the ice point, which corresponds to the values A = 
7r/2 and r] = n/Q. In this case the boundary correlator H^'^ is equivalent to the refined 1- 
enumeration of ASMs. 

In this case the boundary correlator is given by formula (15.101) where 

, , sine , , sinfe — 7r/3) 

sm(£ — 7r/3) sm(£: + 7r/3j 

andp7v-i(a^) is Continuous Hahn polynomial, see (14.171) . The finite difference equation (15.51) 
reads 



1 ixA /2 ix 
3 ~ "6 / is " 6" 



Pn~i{x + 6i) + 



r 4 

— ---NiN -I) 
18 9 ^ ' 



Pn~i[x) 



1 ix\ ( 2 ix\ 
+ o + VN^x{x - 6i) = (5.31) 



3 6 y V3 6 

Similarly to the previous case, employing pAr_i((9e ± 6i) = e^^'^pAr_i((9e)e^^'^ we obtain 



where is the second order differential operator 



= (5.32) 

e=0 



Ke = \ ^\^'iedl sin3£ + isin2 3£-iV(iV-l). (5.33) 
9 9 

Now our aim is to obtain the recurrence relation for and find its solution. 
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The operator in this case reads 



Ds =: ujd^ = -1= sine sin(£: — 7r/3)9e. 

V3 

Taking into account the identity 

sinSe = — 4 sine sin(£: + tt/3) sm{e — n/3) 



and 



sine 



UJ 



2 V^2_^ + i' 



sin(e + 7r/3) 



V3 uj-l 



2 ^oj^-u + l' 



(5.34) 



(5.35) 



(5.36) 



we reexpress the operator in terms of uj and the operator D^, with the latter acting from the 
very left, 



Dl{uj - if - Dsiuj^ - 1) - A^(A^ -l){uj^-uj + 1] 



1 



Taking into account that 



9-\Deg) 



u-l' 



for the operator := g ^^^K^g'^ ^ we obtain 



K, = Di{uj - 1) - D,[{2N - l)uj + 1] + N{N - 1) 



u - 1 



- u; + 1 



(cu^ — c<j + l)/(co' — l)in (15.321) allows us to write 



Choosing T{e) 

PN-i{de) g''-'\Dlico - 1) - D,[i2N - l)co + 1] + iV(iV - 1)^'^-^ 



e=0 



(5.37) 



(5.38) 



(5.39) 



(5.40) 



which, recalling that = ud^, and Eqn. (15.101) . immediately gives us the following recurrence 
relation 



r(r - 2A^+ l)if^ 



(r+l) 



{r-N){N + r- l)H 



(r) 
N 



0. 



This recurrence can be easily solved modulo a normalization constant 



if^^ = const X 



{N + r -2)\{2N -l~r)\ 



(5.41) 



(5.42) 



{r-l)\{N~r)\ 

A possible way to fit the normalization condition (12.101) . is to consider the generating func- 



tion Hn{z) defined via Eqn. (15.151) . The result reads 

_ (2iV-l)!(2iV-2)! 
' (iV-l)!(3iV-2)! 



1-N, N 
2-2N 



(5.43) 
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where the proper normalization is easily determined through Chu-Vandermonde identity 

2fJ~'^' ^ l] =^^Z3^. {a)^:=a{a + l)---{a + m-l). (5.44) 
Inspecting the coefficient of z^'^^ in (15.431) we finally obtain 

r(r) _ I Af -1 / I N -1 J 
(3N-2\ 

The refined 1-enumeration of ASMs, Eqn. (13. 5I) . immediately follows from the last formula and 
the relation (13.81) . 

It is worth to noting that the proof presented here for the refined 1-enumeration of ASMs is 
considerably simpler in comparison to that of Refs. [23,24], which were based on the inhomo- 
geneous square ice partition function formula of Ref. [2]. 

5.4. The A = —1/2 symmetric point and the refined 3-enumeration of ASMs 

(r) 

We shall now apply the same approach to compute the boundary correlator H}^' in the case of 
the A = — 1/2 symmetric point, i.e., when i] = n/3 and A = n/2. This case corresponds to the 
refined 3-enumeration of ASMs. 

As shown in Section|431 in this case there are two sets of Continuous Dual Hahn polynomi- 
als, with differently specified parameters, see Eqns. (14.311) and (14. 321) . which are related to the 
determinant of the Hankel matrix Z, see Eqn. (12.61) . The appearance of two sets of polynomials 
is due to the factorization of the Hankel determinant det Z. The explicit form of this factoriza- 
tion in turn depends on whether is odd or even, see Eqn. (14.261) . Such factorization occurs 
also for the determinant of matrix >V, defined by (15.21) . Denoting D^ix) = det >V, similarly 
to (14.261) we have 

D,,Jx) = Dl^^xDl^\x'), D,^^,ix) = Dt\,{x^)D^^\ (5.46) 

where Dm\x'^) stands for the determinants of the matrices built from the even moments of 
measures with entries of the last column replaced by x"^^ . These determinants are pre- 
cisely the Continuous Dual Hahn polynomials, p~^_i{x'^), specified by Eqns. (14.311) and (14.321) . 
Therefore up to overall constants determined by Dm and , we have 

{const X p^m (x^) if is odd, A^ = 2m + 1 

const X X pm [x ) if A^ is even, A^ = 2m + 2 

In this way the polynomials appearing in (15.31) are expressed in terms of Continuous Dual Hahn 
polynomials. 
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Hence, in the considered case representation (I5.10I) acquires the form 

i7S+2 = const XpW(a2) ^^^2^+1^.-1 



e=0 



H-J,^, = const xpZmg'"''-'-^ 

where 



£=0 



(5.48) 
(5.49) 



, , sine 1 

u; = uj{e) = — — 7^, g = -. (5.50) 

sin(e + 7r/3) to — 1 

Here we have shifted m ^ m + 1 for the odd case for later convenience. 

Starting from expressions (15.481) and (15.491) . one can derive the recurrence relation for H^^^ 

in the cases of even and odd, respectively. Here instead we shall proceed differently using the 

fact that it is possible to express the boundary correlator in terms of a single set of polynomials, 

and thus treat both cases in a unified and simplified way. Indeed, denoting 

u^2li^):=pl:\x'), a = 0,1, (5.51) 

let us consider the polynomials 

/x^ 1 1 1\ 

^-m{x) = Sm[^; j. (5.52) 

These polynomials arise when one studies the action of the 'forward shift operator' for the 
Continuous Dual Hahn polynomials (see, e.g., Ref. [30], Eqn. (1.3.7)), 

U2m+2{x + 3i) - U2m+2{x - 31) = -i(m + 1)| 4^j,(x). (5.53) 

Surprisingly enough, changing the sign in LHS of this relation gives us again Continuous Dual 
Hahn polynomials, which are exactly those defined above as (^)' 

U2m{x + 3i) + U2m{x - 31) = 2ufl{x) . (5.54) 

This relation can be easily proven by expressing all polynomials on both sides as truncated 
hypergeometric series. It is worth to note that this relation is not a specialization of some 
general relation for the Continuous Dual Hahn polynomials, but is instead specific for the 
particular choice of parameters of the polynomials. 

As a direct consequence of relations (15.531) and (15.541) we can rewrite the correlator as 



^^2^+2 = const X U2m+2m sm3e g^^+W'^ 

H2I+3 = const X U2m+2{de) COS 3^ ^2™+ V^^ 



(5.55) 

e=0 

(5.56) 

e=0 
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Noticing that 



sin 3e 



3^3 uj{uj + l) 



cos 3e 



(cu- l)(2t<; + l)(cj + 2) 



(5.57) 



2 (CU2 + CU + 1)3/2' 2(CU2 + + 1)3/2 

and recalling that g = 1/ {lu — 1) it can be easily seen that the correlator possesses the structure 



H. 



(r) 



B. 



(r-l) 
2m 



5, 



(r-2) 



2m 



2m+2 



H. 



(r) 

2m+3 



25, 



2m 



+ 55j„-') + 2B 



(r-2) 
2m 



(5.58) 
(5.59) 



(r) 

where the quantities -Bg^ ^re defined as 

-^2m i>2mU2m+2{dt 



e=0 



{UJ^ + UJ + 1)3/2 

Here, 62m is some normalization constant; we assume that 

2m 

E^22 = i. 



r = 0,1, 



2m. 



(5.60) 



(5.61) 



r=0 



This condition, together with Eqns. (15.581) and (15.591) . ensures that normalization condition 
(12.101) is satisfied. As in previous cases, the proper normalization will be restored at the end of 
computation, according to condition (15.611) . 

Thus, instead of studying the correlator for N even and odd separately it is enough to 
consider the quantity i?2m which is defined by an essentially similar formula, Eqn. (15.601) . The 

(r) 

procedure developed previously will be applied now to derive a recurrence relation for i?2m- 
The finite-difference equation (15.71) reads 



(1 + X^) M2m+2(3: + 61) - U2m+2{x - 61) - 24i(m + l)xU2m+2{x) 

Using this equation, like in the previous cases, we can write 



0. 



U2m+2ide)KJ{e] 







e=0 



(5.62) 



(5.63) 



sin Be cos 3e (9^ + l) — 6(m + l)^^ 



(5.64) 



where is the second order differential operator 

while f{e) is some arbitrary function. Obviously, to obtain a recurrence relation for the quan- 
tities -B2m' the function f{e) is to be chosen of the form 



f = r 



g2m+l^r+2 
(1 + + Cj2)3/2 
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(5.65) 



where again r is some function to be chosen later. 

Taking into account that the differential operator in the considered case reads 

2 



and passing to the variable u, we obtain 



sine sin(£: + Tr/S)^^ 



K, = D 



2(2cj2 + 5w + 2)(cj2- 1) 



UJ{UJ^ + UJ+1) 

+ 4w + l)(2cj2 + + 2)(cj - If 



ujiuj^ + uj + lV 



+ Am 



(Cj2 + CJ + 1)' 



- + Am 



4u;(cj2 + + 1) 



UJ 



Defining operator by the formula 

-,2m+l 



9 



9 



2m+l 



{UJ^+UJ + 1)3/2 (cu2 + + 1)3/2 

we may rewrite condition (15.631) 

-,2m+l 



[UJ^ + UJ + 1)'^'^ 



r+2 



e=0 



where 



(5.66) 



(5.67) 



(5.68) 



(5.69) 



= i (2u;^ + 5u + 2){lu^ - 1) 



{2u? + 5cj + 2)(7c<;2 - 2cj - 1) + 4m(cu^ + 5cu^ + 4c<;2 - 1) 



+ 2(2a;^ + 5a; + 2)(5a;^ - 2a; + 1) + 2m(4cj^ + 260;^ + 17a;^ + 5a; + 2) 



+ Am^u){?iU? + 4a; + 2) 



1 



a;(a;2 + a; + 1) 



. (5.70) 



Choosing r = a; + l + a; ^we therefore obtain 

,2m+l 



U2m+2{de 



9 



2 + 1)3/2 



a; 



D'^ (2cu^ + 5cj3 - 5cj - 2) 



(14tu^ + Slu^ + 2uj^ -9uj-2)+ 4m(cj^ + 5uj^ + Au? - 1) 



+ 2(10cj^ + 21cj^ + 2cj2 ^ + 2) + 2m(4cj^ + 26cu^ + YIu? + 5cj + 2) 

+ 4m2(3cj3 + 4cj2 + 2cu)lcu' 

J £=0 



0. (5.71) 
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(r) 

This equation leads immediately to the recurrence relation for -B2m- The latter enjoy the sym- 
metry i?2m = B^m~^^ which obviously guarantees H^'' = H^~^'^^^ (this is also known as 
the top-bottom or left-right symmetry of the set of ASMs). To make this, symmetry more ap- 
parent it is convenient to write this recurrence relation in terms of Em ■= B^~^^\ so that 

(r) (—r) 

r = — m, . . . , m. The recurrence relation reads 

2(r - m - 2)(r + m + 1) E^^~^^ + (Sr^ + lOrm + r - Sm^ - 9m - 6) 
+ 2(1 + 8m)r - (5r^ - lOrm - r - Sm^ - 9m - 6) 

- 2(r - m - 1) (r + m + 2) E^^^^ = 0. (5.72) 

(r) (r) 

Using this relation one can find recursively all Em s for any given m assuming that Em = 
if |r| > m. However, since the recurrence relation is five-term it can hardly be solved, e.g., 
by guessing its solution. The remaining of this Section is an exposition of a possible way to 
solve it explicitly. As we shall show now this can be done by successive transformations of the 
generating function for B2m. 



Consider the generating function 



2m 



Em{z) := Yl = ^'™E^2nl^^ Em{z) = Em{z-'). (5.73) 

r=—m r=0 

which satisfies, as a direct consequence of recurrence relation (15.721) . the following homoge- 
neous second order linear differential equation: 



[z - z-^){2z + 5 + 2z-^){zd,f 



zd,. 



(22 + 5 + 2z-^){'iz - 2 + 3^-^) + 2m(5;z + 8 + ^z-^] 

m{6z - 1 + 6z-^) + m^{2z + 3 + 2^"^) = . (5.74) 



-(z- z-^ 



The solution of this equation we are interested in is a polynomial of degree 2m times factor 
2;"'", see Eqn. (15.731) . Let us now consider the substitution 

^ _ — ^ _ gxp(i7r/3) (5.75) 
qx — 1 

which maps the six singularities of this equation lying on the real axis of the complex 2; -plane, at 
points z = —1, 0, 1/2, 1, 2, 00, onto the six roots of equation x^ = 1, i.e., x = l,q,q'^, q^, g^, g^, 
lying on the unit circle of the complex a; -plane. Note moreover the trivial but useful identity 
1 + q^ = q. Taking into account symmetry (15.731) enjoyed by Em{z), and the identity {qx — 
l)(x — q) = {qx){x — 1 + x~^), it is easily seen that the function 

Vm{x) := (x - 1 + x-^r Em f-^^l (5.76) 

\ qx ~ I J 
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is again of the form x^™ times a polynomial of order 2m in x, symmetric under x 1/x, i.e., 
the function V^m(e"^) is an even trigonometric polynomial of degree m in ip. Equation (15 .741) 
translates into the following equation: 



+ 



(x + 1 + x"^)(x^ - 5x + 12 - 5x"^ + x"^) - 2m(x + x~^)(x^ - 5 + x~^) 



Xdr 



— X — X 



m(x2-4x + 13-4x"^+x~^)-m^(x^-7 + x~2) lv;„(x) = 0. (5.77) 



At the first glance there is no advantage in considering the function Kn(x) since the underlying 
recurrence relation for expansion coefficients of Vm{x) is even worse than that for those of 
Ern{z): it is a seven-term relation. However, if one considers instead the function. 



/1„(X) = (X - X-^)2"^+^ (X + 2 + X-')Vrn{x) 



(5.78) 



then it appears that the differential equation satisfied by h„i{x) contains only integer powers of 
x^ in their coefficients 



x-'-x'''){xd^y~ 3(2m+l)(x''+x-'')-6 x9^ + (3m + l)(3m+2)(x''-x"'') }hmix) = 



(5.79) 



or, using the variable 0, related to x as x = exp{i(p), we have 

1 



(2m + 1) cot 3(/9 — 



- (3m + l)(3m + 2) )h^{e'^) = 0. (5.80) 



sin3(/? 

It is easy to see that the underlying recurrence relation now is just a three-term and it appears 
to be solvable explicitly. 

The particular solution of Eqn. (15.801) which we are looking for is the function h„t{e^'^) being 
an antisymmetric trigonometric polynomial of degree 3m + 2 in see Eqn. (15.781) . Hence, we 
are led to use the substitution of the form 



2m+l 



/i(e'^) = ^ 7fesin[(3m + 2 - 3A;)y?]. 



(5.81) 



fe=0 



The coefficients 7^ have to satisfy the recurrence relation 

(3fc + 2)(3A; + 3)7fc+i + (3m + 2-3A;)7fc- [3(2m-A;) + 7] [3(2m- A;) + 6]7fc_i = 0. (5.82) 

Inspecting the explicit form of the first few 7fc's, allows us to guess that the solution of the 
recurrence relation is 



72« = 7o (-1) 



^ (-m-2/3)i/m^ i {-m - 2/3)i+i f m 

I21+1 = 7o (-Ij 



(1/3). 



(1/3) 



(5.83) 



27 



that can be easily verified directly. Here {a)i stands for Pochhammer symbol, defined previ- 
ously in Eqn. (15.441) . Thus, we have obtained for the function hm{x) the following expression: 



hm{x) 



where the functions fmix) and gmix) are given by 

2/3 



3m + 2 



3m + 1 



ix. 



(5.84) 



9m[X) 



fm{x) 



E 

fc=0 
m 

E 

k=0 



m 



k 



^ ^/"^^ ^^3m+2-6fc _ ^-3m-2+6fc 



m + 1/3 



m — k 

m — 1/3 
m — k 



) 



X' 



3m+l— 6A: 



— X 



-3m— 1+6A: 



) 



(5.85) 
(5.86) 



and Cm is some constant such that V"m(l) = Em{l) = 1; this normalization follows from the 
normalization condition (15.611) for -B2m- 

It is worth noticing that function hm{x) as given in Eqns. (15. 841) . (15.851) . (15.861) . in connec- 
tion with the problem of refined 3-enumeration has also been found by Stroganov in Ref. [27], 
within a different approach, using a certain functional equation satisfied by the inhomogeneous 
square ice partition function. This functional equation had already been investigated in [31] as 
the Baxter T-Q equation for the ground state of XXZ Heisenberg spin-1/2 chain at A = —1/2 
and with odd number of sites N = 2m + 1. 

The problem we are facing now is to reconstruct, from the explicit knowledge of function 
hm{x), the generating function Em{t). This amounts essentially to divide out the factors (x — 
^-i-j2m-i-i ^j^^ (^3; _|_ 2 -|- x~^) to find function Vm{x) first, and next to change back to the original 
variable z = (x — q)/{l — qx) to recover Emiz). We shall follow the line of our previous 
paper [32] where this procedure was fulfilled. 

To undertake the first step in this program it is useful to note that the functions fm{x) and 
gm{x) can also be written as follows 



fm{x) 



r(m + 1/3) 
m!r(l/3) 



r(m + 4/3) 
m!r(4/3) 



X 



3m+2 



-m . 



X 



-m - 2/3 
1/3 



-3m— 2 TP 
2-fl 



X 



-m . 



-m - 2/3 
1/3 



X 



(5.87) 



X 



-3m-|-l 



1/3 



4/3 



X 



X 



3m— 1 



'Fi 



-m . 



-m + 1/3 
4/3 



X 



(5.88) 



Since the parameters of the hypergeometric functions entering these expressions differ by in- 
tegers one can expect that gm{x) and fm{x) are connected by some three-term relations via 
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Gauss relations (see, e.g., §2.8 of [29]). Indeed, using Gauss relations it can be shown that 

3m + 4 



2^1 



-m , — m — 2/3 
1/3 



m — 1 , —m — 2/3 
4/3 

3m + 2 , ^ 
(1 + C) 2F1 



m , — m + 1/3 
4/3 



and therefore one can express the function gmix) in terms of fmix) and fm+i{x 

9mix) 



3m + 2 n „ , - 3(m + 1) 

(X^ + X ^) - ^ , fm+l[X). 



2(3m+ 1) ' ' ' 2(3m+ 1) 

Substituting (15.901) into (15.841) we obtain an analogous formula for function hm{x) 



C (5.89) 



(5.90) 



3m + 2 
2(3m + 1) 



3m + 2 



Introducing now the function Qm{^), implicitly defined by 



fm{x) 



[X — X 



-l\2m+l 



(5.91) 



(5.92) 



the factor (x — x i)2™+i (x + 2 + x ^) can be formally extracted in expression ( 15.911 ) for h„i{x), 
thus giving us a representation for Vmix) in terms of Qm{x) and Qm+i{x), 



^m{x^ '-'1 



3m + 2 
2(3m + 1) 



_|_ 2 

(x - 1 + x"^)^ (5m(a;) - — (x - 2 + x~^) Qm+i{x) 

im + 2 



(5.93) 

The meaning of this procedure becomes apparent by noticing that the function Qm (x) can 
be found explicitly from expression (15.881) for the function fm{x) in virtue of the so-called 
cubic transformation for the Gauss hypergeometric function. The details of this derivation are 
given in Appendix A. For function Qm{x) the following explicit formula is valid 

(2m)! / qx^^ — q^^x^ ™ 



Qm (3') 



3™ {m\y 



q-q 



-1 



-m , m + 1 
—2m 



qx — q ^x ^ 



gx"^ — q ^x 



(5.94) 



Hence, function Vm{x) can be found by inserting this expression into (15.931) . thus giving 

(2m)! (2m + 1)! 



Vmix) 



ml (3m + 1)! 



fx — 1 + X ^ 



-1 -1 \ ^ 
^2 I QX — q X ^ 



-m , m + 1 
—2m 



qx — q -^x ^ 



qx'^ — q ^x 



2(2m + 1) 
3m + 2 



(x — 2 + x^ 



— 1 —1 \ rn+l 

/ gx — g X ^ 



X 2F1 



q-q 
m — 1 , m + 2 
-2m - 2 



qx — q ^x ^ 
qx~^ — q~^x 



(5.95) 
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Here we have written the expression for Vm{x) taking into account also the proper normal- 
ization of this function, = 1. The normalization can be verified by virtue of Chu- 
Vandermonde identity (15 .441) and it is equivalent to the choice 



(3m + i; 



gm+l^l (2m + 2)! 
(3m + 3)! 



(5.96) 



of this constant in Eqns. (15.841) and (15.931 ). It is worth to mention that functions Qmi^), 
Eqn. (I5.94I) . and Vmix), Eqn. (15.951) . are the first and the second solution, respectively, of 
Baxter T-Q equation for the ground state of XXZ Heisenberg A = —1/2 spin chain with odd 
number of sites N = 2m + 1, see Refs. [24,27,31,33]. 

To obtain function Em{z) from the given expression for Vmix) one can use the formula 



(z + l + Z-'rVrr, 



z + q 
zq + 1 



(5.97) 



which is the inverse of (15.761) . Applying this transformation to (15.951 ) we obtain 



Er, 



_ (2m)! (2m + 1)! 1 
~ 3'"m! (3m+ 1)! (z + 1 + z-^)'^ 
2(2m + 1) , 



9(j + 2)'"2Fi 



m+1 



-m , m + 1 
—2m 
m — 1 , m + 2 
-2m - 2 



z 



-1 



+ 2 



z 



z + 2 
1 + 2- 



z + 2 



(5.98) 



This expression is however not the final answer yet, since the factor (z + 1 + z~^Y in the 
denominator is to be cancelled explicitly (recall that Em{z) is to be z^'"^ times a polynomial 
of degree 2m in z, hence the expression in the brackets contains implicitly the proper factor 
{z + l + z^^Y). This can be achieved using again Gauss relations; for details, see Appendix B. 
The final result is 



E„ 



(2m)! (2m + 2)! 
3™ (m+1)! (3m + 2)! 



(2m + l)(z + 2)"^ 2Fi 



3m iz + 2)'"-! 2^1 



m , m + 2 
-2m - 1 
m + 1 , m + 2 
—2m 



1 + 2^ 



z{z + 2) 
l + 2z 



z{z + 2) 



(5.99) 



This completes the derivation of function Em{x). 

(r) 

We can now extract from this expression closed formulae for the coefficients -B2m' which 
are related to Em{z) through Eqn. (15.731) . To find these coefficients we expand (15.991) in power 
series in z, thus expressing E^iz), as a triple sum, and next we apply Chu-Vandermonde 
formula (15.441) to make the sum with respect to the index defining the hypergeometric series 
in (15.991) . thus expressing E.m{z) as a double sum. These two summations can be rearranged 
in such a way that one of them becomes with respect to r while the other one defines the 



30 



coefficients of power expansion in z. We obtain 



B. 



(r) 
2m 



(2m + l)!m! 
3*^ (3m + 2)! 



[r/2] 



(2m + 2 - r + 2£) 



/3m + 3 



<?=max(0, r—m) 



X 



2m + i-r + 1 
m + 1 



m 



m + 1 



ir-2f 



(5.100) 



Here [r/2] denotes integer part of r/2. This expression for indeed solves the five terms 



,(m+r) 
2m 



recurrence relation (15.721) (recall that Em = B 

We would like also to mention that formula (15.1001) can also be written in terms of termi- 
nating hypergeometric series, for instance, as follows: 



B. 



(r) 
2m 



2r /3m+3\ /2m+l-r\ 
V r J \ m+1 / 

{ m+1 ) 



2 4^3 



-(r- l)/2, -r/2, m + 2, 2m + 2-r 
(3m + 4 — r)/2 , (3m + 5 — r)/2, m — r + 1 



m + 1 



4-^3 



-(r-l)/2, -r/2 + 1, m + 2,2m + 2-r 
(3m + 4 — r)/2 , (3m + 5 — r)/2, m — r+1 



(5.101) 



2m 



This formula is valid for r = 0,1,. ..,m (a similar expression forr = m + l,m + 2, 
can be simply obtained through the replacement r 2m — r in RHS of (15.1011) . The two 4F3 
in (15.1011) can be further combined into a single 5F4 . Analogous formulae for i?2m in terms 
of terminating hypergeometric series of argument 4 may be written down as well. Analyzing 
these expressions, however it seems to be hard to perform the sum in (15.1001) in a closed form, 
even if very suggestive similarities can be found with known summation formulae, see §§7.5 
and 7.6, especially §7.6.4, of Ref. [34]. 

(r) 

The complete expression for the boundary correlator Hj^ can be readily obtained by in- 
serting (15.1001) or (15.1011) into (15.581) and (15.591) . Finally, the refined 3-enumeration of ASMs, 
A(N, r; 3), follows from multiplying the result for H^^ by the total number of 3-enumerated 



ASMs, A{N, 3), see (ESJ) and (|33 



6. Conclusion 

The main purpose of the present paper was to point out the close relation between ASM enumer- 
ations and some classical orthogonal polynomials. In particular, in Section|4|we have shown an 
alternative way to recover known results for the partition function of the six- vertex model with 
DWBC, in the three cases of the A = line, and the symmetric A = 1/2 and A = — 1/2 points, 
corresponding to 2-, 1- and 3-enumeration of ASMs, respectively. The derivation we have pre- 
sented in Section m is in our opinion extremely simple and straightforward. It is based on the 
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fact that the Hankel determinant entering the representation for the partition function can be 
naturally related in these three cases to Meixner-PoUaczek, Continuous Hahn, and Continuous 
Dual Hahn polynomials, respectively. 

It is to emphasize that the three considered cases, x = 1,2,3, are the only ones in which 
'factorized' answers for ASMs' x-enumerations are known to exist. The fact that no set of poly- 
nomials corresponding to other choices of the parameters can be found in the Askey scheme 
strongly suggests that no 'factorized' answer exist for other values of x. However, one might 
speculate that for some values of the crossing parameter r], the partition function of the six- 
vertex model with DWBC might still admit a 'factorized' form, but in terms of g-numbers 
(with q related to t] through q = e^"'), and obtained via some suitable g-polynomial, possibly of 
Askey-Wilson type. Concerning this, it is worth mentioning a recent paper [35] which contains 
rather promising preliminary results in this direction. Speculating further, an interesting prob- 
lem which we would like just to hint here concerns a possible relation of such g-polynomials 
as q tends to a cubic, quartic or sixth root of unity, with the three sets of classical polynomials 
mentioned above. In this respect, recall that in paper [23] g-Legendre polynomials were used 
while in the present paper Continuous Hahn polynomials have shown to play an analogous role. 

In Section |5] we have used the knowledge of appropriate orthogonal polynomials for the 
cases under investigation to derive recurrence relations for the boundary one point correlator 
H^^\ In the first two cases, such recurrences are trivially solved, and known result for the 
closely related problems of ASMs'refined 1- and 2-enumerations are easily reproduced. The 
same approach is applied to the symmetric A = —1/2 point, but the resulting recurrence rela- 

(r) 

tion for H}^' appears very intricate, not being two-term (like in cases of 2- and 1-enumerations) 
but rather five-term. Even though the differential equation for the corresponding generating 
function is not of hypergeometric type, it has been shown to be solvable in terms of a suit- 
able linear combination of hypergeometric functions. The so-called cubic transformation for 
the Gauss hypergeometric function has been applied to work out an explicit formula for the 
refined 3-enumeration of ASMs, A{N, r; 3). The latter appears not to be writable as a single 
hypergeometric term, i.e., not to be 'round' (or 'smooth'), contrarily to other known expres- 
sions for enumerations of ASMs, and it hardly could have been conjectured on the basis, e.g., 
of computer experiments. 
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Appendix A 



We reproduce here for the sake of completeness the proof of Eqn. (15.941 ) given in a previous 
paper [32]. 

The key identity which is to be used here is the so-called cubic transformation of Gauss 
hypergeometric function [29] which in its most symmetric form reads: 



r(a) /a + 1/3 , a 

2-^1 



r(2/3) 



2/3 



3> _i r(a + 2/3) /a + l/3,a + 2/3 
^ '"'^ ^ r(4/3) 4/3 



-3a+l 



1-C 

1-UJ 



—3a 



r(3a) 



2Fl 



a + 1/3 , 3a 



UJ- 



1-C 



. (A.l) 



r(2a + 2/3) V 2a + 2/3 

Here to is a primitive cubic root of unity, to = exp(±2i7r/3), and a is arbitrary parameter. To 
show that indeed the cubic transformation is relevant to our case, let us rewrite (15.881) in the 
form consistent with LHS of (lA.ll) . Taking into account that 



2^1 



and 



-m , — m + 1/3 



4/3 



r(l/3)r(4/3) 



r(-m + l/3)r(m + 4/3) 



-m , — m — 1/3 
2/3 



r(l/3) 



-ir(-m-l/3) 



r(m + 4/3) 



r(m + 4/3) ' ' r(2/3) ' r(-m + l/3) 

it is easy to see that (15.881) can be rewritten in the form 

;-l)™+i(3m + 1) 



-1^ 



33m+l ^! 



(A.2) 
(A.3) 



fm{x) 



33m+l ^^n2 



X 



3m+l 



+ X~ 



r(— m — 1/3) f —m , —m — 1/3 

f(273) ' ' 
^r(-m + l/3) 



2/3 

-m , ~m +1/3 
4/3 



X 



X 



(A.4) 



r(4/3) 

Clearly, both terms in the brackets are the same as in LHS of (lA.ll) provided the parameter a is 
specialized to the value a = —m — 1/3. 

To apply the cubic transformation to (IA.4I) we first define 



W{a-Q :-- 



so that 



r(a) ^a + 1/3 , a 



r(2/3) 



2/3 



3, r(a + 2/3) /a + l/3,a + 2/3 
^ ' ^ r(4/3) ' 'V 4/3 



'-1 



\m+l 



(3m + 1)! 



33m+l (^1)2 



X 



3m+l 



W{-m - 1/3; x-^). 



(A.5) 
(A.6) 
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Next, we note that for a sum of two terms one can always write 



X + Y 



-1 



(A.7) 



q-q - ^ q-q- 

and if g = exp(i7r/3), which is exactly the case, one can set uj = for the first pair of terms 
and oj = for the second one. This recipe allows one to apply the cubic transformation, that 
gives 



W{a-0 



-3a+l 



r(3a) 



C) 



-3a 



r(2a + 2/3) g(l-g2)-3a+i 



2^1 



a + 1/3 , 3a 
2a + 2/3 



I ^3a+l 17 
+ q 2-ri 



1-C 

a + 1/3 , 3a 
2a + 2/3 



' X-q~' 
i-C 



. (A.8) 



To obtain a new formula for fm{x) via (IA.6I) we have to evaluate now the limit a — m — 1/3 
of (IA.8I) . The limit of the pre-factor can be easily found due to 

r(3a) 2 (-l)'"+i(2m)! 



a_l™ 1/3 r(2a + 2/3) 3 (3m + 1)! 
To find the limit of the expression in the brackets in (IA.8I) we note that 



2C-q' 
q r =: s 



1-C 



q — — = I — s 



1-C 



(A.9) 



(A. 10) 



and hence the following formula can be used 

2F1 



lim 

a— »— m— 1/3 



a + 1/3 , 3a 
2a + 2/3 



— - 2-ri 



-m , —3m 
—2m 



s). (A.ll) 



Formula (lA.l II) can be proved, for instance, by virtue of standard analytic continuation formu- 
lae for the hypergeometric function (see, e.g., Eqns. (1) and (2) in §2.10 of [29]). Collecting 
formulae we arrive to the expression 



W{-m -1/3;C) 



(2m) 



C) 



3m+l 



2^1 



m , —3m — 1 
—2m 



(3m + 1)! (q — q ^] 
Finally, substituting this expression into (IA.6I) and using the identity 



1-C 



a , b 



C =(i-C)-%i^i 



a , c — b 

c 



c 



c 



(A.12) 



(A.13) 



we obtain 



(2m) 



3™(m! 



[X — X 



-l\2m+l 



qx — q x^ 



q-q 



-1 



-m , m + 1 
—2m 



1^-1 



qx — q X 



qx ^ — q ^x 



(A. 14) 



Obviously, this expression leads directly to (15.941) which is thus proved. 
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Appendix B 

Here we explain how the final expression for function Em{z), Eqn. (I5.99I) . can be obtained 
from Eqn. K9^ . 

We begin with noting that the hypergeometric functions (polynomials) which enters expres- 
sion (15.981 ) belong to the class of the functions 

a , b 
^a-b + 1 

which are, in the case of a being a negative integer, possess the symmetry 



(B.l) 



2F1 



a , b 
a-b+1 



a , b 
a-b+1 



(B.2) 



Obviously, since Em{z) is symmetric with respect to z — > 2;"^, it is natural to deal only with 
the functions of the form (IB. II) when transforming the expression for Em{z). It is convenient 
to use the notation 



u) 



-m , k + 1 
—m — k 



z-^ + 2 
z + 2 



u:= Z + 1 + z~^. 



(B.3) 



Note that '^^m{u) is a polynomial of degree m in u. Using Gauss relations the following 
identities can be proven 

m{m + 2/c + 1) 



m \ 



U] 



m + 2k + 2 



{m + k + l){m + k) 



m 



(2n + 3)v[/«^(^) (B.4) 
'u + ?,)^!tt'Xn) . (B.5) 



2(m + A; + 1) * ™ ' ' ' 2{m + k + 1] 
Returning to (15.981) we note that in terms of ^m^'s the generating function Em{z) reads 



(2m)!(2m+l)! 1 
(3m + 1)! 

Applying (IB. 41) . with k = m + 1, gives us 



Emiz) 



(2m)!(2m+l)! f 2(2m+l) ^(^^^^ ^ 



3'"m! (3m + 1)! [ 3m + 2 



6m 



(m+l) 
m— 1 ' 



U 



+ 



9 



m 



. 2(2m + l) 
3m + 2 ™ ^ ^ 3m + 2 



3m + 2 



u 



■ (B.7) 



Now, the relation (IB.5I) with k = m shows that the expression in the brackets is zero. Hence 
we find 

(2m)! (2m + 2)! 



Err,(z) 



i2m+l)¥::^-^'\u)-3m¥Zriu) 



3™ (m + l)! (3m + 2)! 
Restoring the original notations, we arrive to expression (15.991 ). 



(B.8) 



35 



References 



[1] V.E. Korepin. Calculations of norms of Bethe wave functions. Commun. Math. Phys. 
(1982) 86 391^18. 

[2] A.G. Izergin. Partition function of the six-vertex model in the finite volume. Sov. Phys. 
Dokl. 32(1987) 878-879. 

[3] A.G. Izergin, D.A. Coker, and V.E. Korepin. Determinant formula for the six-vertex 
model. /. Phys. A: Math. Gen. 25 (1992) 4315^334. 

[4] N.M. Bogoliubov, A.G. Pronko, and M.B. Zvonarev. Boundary correlation functions of 
the six-vertex model. J. Phys. A: Math. Gen. 35 (2002) 5525-5541. 

[5] V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin. Quantum Inverse Scattering Method 
and Correlation Functions. Cambridge University Press, Cambridge, 1993. 

[6] W.H. Mills, D.P. Robbins, and H. Rumsey. Proof of the Macdonald conjecture. Invent. 
Math. 66 (1982) 73-87. 

[7] W.H. Mills, D.P. Robbins, and H. Rumsey. Alternating-sign matrices and descending 
plane partitions. J. Combin. Theory Sen A 34 (1983) 340-359. 

[8] D.P. Robbins and H. Rumsey. Determinants and alternating-sign matrices. Advances in 
Math. 62 (1986) 169-184. 

[9] D. Zeilberger. Proof of the alternating sign matrix conjecture. Elec. J. Comb. 3 (2) (1996) 
R13. 

[10] G. Kuperberg. Another proof of the alternative-sign matrix conjecture. Internat. Math. 
Res. Notices 1996 (1996) 139-150. 

[1 1] D.M. Bressoud. Proofs and Confirmations: The Story of the Alternating Sign Matrix Con- 
jecture. Cambridge University Press, Cambridge, 1999. 

[12] N. Elkies, G. Kuperberg, M. Larsen, and J. Propp. Alternating-sign matrices and domino 
tilings (part I). /. Algebraic Combin. 1 (1992) 111-132. 

[13] N. Elkies, G. Kuperberg, M. Larsen, and J. Propp. Alternating- sign matrices and domino 
tilings (part II). J. Algebraic Combin. 1 (1992) 219-234. 

[14] H. Cohn, N. Elkies, and J. Propp. Local statistics for random domino tilings of the aztec 
diamond. Duke Math. J. 85 (1996) 117-166. 



36 



[15] W. Jockush, J. Propp, and P. Shor. Random domino tilings and the arctic circle theorem. 
Preprint, arXiv: .math.CQ/9801068 

[16] A.V. Razumov and Yu.G. Stroganov. Spin chains and combinatorics. 7. Phys. A: Math. 
Gen. 34 (2001) 3185-3190. 

[17] P. Di Francesco and P. Zinn- Justin. Around the Razumov-Stroganov conjecture: proof of 
a multi-parameter sum rule. Preprint, arXiv: math -ph/04 1 006 1 

[18] E.H. Lieb. Exact solution of the F model of an antiferroelectric. Phys. Rev. Lett. 18 (1967) 
1046-1048. 

[19] E.H. Lieb. Exact solution of the two-dimensional Slater KDP model of a ferroelectric. 
Phys. Rev. Lett. 19 (1967) 108-110. 

[20] B. Sutherland. Exact solution of a two-dimensional model for hydrogen-bonded crystals. 
Phys. Rev Lett. 19 (1967) 103-104. 

[21] E.H. Lieb and FY. Wu. In: Phase Transitions and Critical Phenomena, Vol. 1, edited by 
C. Domb and M. S. Green, Academic Press, London, 1972, p. 321. 

[22] R.J. Baxter. Exactly Solved Models in Statistical Mechanics . Academic press, San Diego, 
1982. 

[23] D. Zeilberger. Proof of the refined alternating sign matrix conjecture. New York J. Math 2 
(1996) 59-68. 

[24] Yu. G. Stroganov. A new way to deal with Izergin-Korepin determinant at root of unity. 
Preprint, arXiv: math-p h/0204042, 

[25] V.E. Korepin and P. Zinn-Justin. Thermodynamic limit of the six-vertex model with do- 
main wall boundary conditions. J. Phys. A: Math. Gen. 33 (2000) 7053-7066. 

[26] P. Zinn-Justin. Six-vertex model with domain wall boundary conditions and one-matrix 
model. Phys. Rev. E 62 (2000) 3411-3418. 

[27] Yu. G. Stroganov. 3-enumerated alternating sign matrices. Preprint, arXiv: 
|:math-ph/030400i| 

[28] G. Szego. Orthogonal Polynomials. Fourth edition. Colloquium Publications, Vol. 23, 
Amer. Math. Soc, Providence, RI, 1975. 



37 



[29] A. Erdelyi. Higher Transcendental Functions, Vol. I. Robert E. Krieger Publishing Com- 
pany, Malabar FL, 1981. 



[30] R. Koekoek and R.F. Swarttouw. The Askey-scheme of hypergeometric orthogonal poly- 
nomials and its g-analoque. Report no. 98-17, Delft University of Technology, 1998. 

[31] Yu. G. Stroganov. The importance of being odd. J. Phys. A: Math. Gen. 34 (2001) L179- 
L185. 

[32] F. Colomo and A.G. Pronko. On the refined 3-enumeration of alternating sign matrices. 
Preprint, arXiv: math-ph/0404045 

[33] G. P. Pronko and Yu. G. Stroganov. Bethe equations 'on the wrong side of equator'. 
/. Phys. A: Math. Gen. 32 (1999) 2333-2340. 

[34] A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev. Integrals and series, Vol. 3: More 
special functions. Gordon and Breach Science Publishers, New York, NY, 1990. 

[35] M. E. H. Ismail, S. Lin, and S. Roan. Bethe anstaz equations of XXZ model and g-Sturm- 
Liouville problems. Preprint, arXiv: math-ph/0407033 , 



38 



